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Abstract 

The paper addresses the problem of the transverse force (Mag- 
nus force) on a vortex in a Gahlean invariant quantum Bose hquid. 
Interaction of quasiparticles (phonons) with a vortex produces an ad- 
ditional transverse force (lordanskii force). The lordanskii force is 
related to the acoustic Aharonov-Bohm effect. Connection of the ef- 
fective Magnus force with the Berry phase is also discussed. 

1 Introduction 

In classical hydrodynamics it has long been known that if the vortex moves 
with respect to a liquid there is a force on the vortex normal to the vortex 
velocity 0]. This is the Magnus force, which is a particular case of a force 
on a body immersed into a liquid with a flow circulation around it (the 
Kutta-Joukowski theorem). An example of it is the lift force on an airplane 
wing. 



1 



The key role of the Magnus force in vortex dynamics became clear from 
the very beginning of studying superfluid hydrodynamics. In the pioneer ar- 
ticle on the subject Hall and Vinen defined the superfluid Magnus force as 
a force between a vortex and a superfluid. Therefore it was proportional to 
the superfluid density ps- But in the two-fluid hydrodynamics the superfluid 
Magnus force is not the only force on the vortex transverse to its velocity: 
there was also a transverse force between the vortex and quasiparticles mov- 
ing with respect to the vortex. The transverse force from rotons was found 
by Lifshitz and Pitaevskii ^ from the quasiclassical scattering theory. Later 
lordanskii Q revealed the transverse force from phonons. The analysis done 
in Ref. (see also Refs. P, 0) demonstrated that the Lifshitz-Pitaevskii 
force for rotons and the lordanskii force for phonons originate from interfer- 
ence between quasiparticles which move past the vortex on the left and on 
the right sides with different phase shifts, like in the Aharonov-Bohm effect 
Since the phase shift depends on the circulation which is a topologi- 
cal charge for a vortex, this is a clear indication of connection between the 
transverse quasiparticle force and topology. 

Later on the analogy between wave scattering by vortex and electron scat- 
tering by the magnetic-flux tube (the Aharonov-Bohm effect) was studied 
in classical hydrodynamics for water surface waves (the acoustic Aharonov- 
Bohm effect jl^]). Scattering of the light by a vortex also results in the 
Aharonov-Bohm interference (the optical Aharonov-Bohm effect As 
follows from Ref. P], the Aharonov-Bohm interference always produces a 
transverse force on the vortex, or the fluxon. For the original Aharonov- 
Bohm effect this force is discussed by Shelankov |jl2| . 

The Magnus force on a vortex in a superconductor was introduced by 



Nozieres and Vinen ||T^. The total transverse force on a vortex is respon- 
sible for the Hall effect in the mixed state. In superconductors not only 
quasiparticles, but also impurities produce an additional transverse force on 
the vortex. A reader can flnd discussion of this problem in the review by 
Kopnin p^. 

Despite a lot of work done to understand and calculate the Magnus force, 
it remained to be a controversial issue. Ao and Thouless pointed out a 
connection of the Magnus force with the Berry phase [16| which is the phase 
variation of the quantum-mechanical wave function resulting from transport 
of the vortex round a close loop. From the Berry-phase analysis Ao and 
Thouless concluded that the effective Magnus force is proportional to the 



2 



superfluid density, and there is no transverse force on the vortex from quasi- 



particles and impurities |T5|, |T7|, [T8|. Such conclusion disagreed with the pre- 
vious calculations of the Magnus force in superfluids and superconductors, 
and therefore generated a vivid discussion [|^, [2^, |21| . 

The present paper addresses this problem p2| which is very important for 
many issues in modern condensed matter physics, field theory, and cosmology. 
The analysis is based on studying momentum balance in the area around a 
moving vortex without using any preliminary concept of a force. The word 
force is a label to describe a transfer of momentum between two objects. A 
careful approach is to give these labels to the various terms in the momentum 
balance equation only after derivation of this equation. I restrict myself with 
the problem of the Galilean invariant quantum Bose liquid. For a weakly 
nonideal Bose gas one can use the Gross-Pitaevskii theory |2^. In this 
theory the liquid is described by the nonlinear Schrodinger equation. At 
large scales the nonlinear Schrodinger equation yields the hydrodynamics 
of an ideal inviscous liquid. In presence of an ensemble of sound waves 
(phonons) with the Planck distribution, which is characterized by a locally 
defined normal velocity Vn (the drift velocity of the Planck distribution), 
one obtains two-fluid hydrodynamics. Eventually the problem of the vortex 
motion in presence of the phonon normal fluid is a problem of hydrodynamics. 

The paper starts from discussion of the Magnus force in classical hydrody- 
namics (Sec. 0). In Sec. ^I define the superfluid Magnus force and the force 
from quasiparticles scattered by a vortex. Connection between the Gross- 
Pitaevskii theory and the two-fluid hydrodynamics is discussed in Sec. ^. In 
Sec. I I discuss scattering of a sound wave by a vortex and show that the 
standard scattering theory fails to give a conclusive result on the transverse 
force because of small-scattering-angle divergence of the scattering ampli- 
tude. The next Sec. |] gives a solution of the sound equation around the 
vortex, which is free from the small-angle divergence. Using this solution 
in the momentum balance one obtains the equation of vortex motion, which 
contains the lordanskii force. The momentum transfer responsible for the 
lordanskii force occurs at small scattering angles where a phenomenon anal- 
ogous to the Aharonov-Bohm effect is important: an interference between 
the waves on the left and on the right from the vortex with different phase 
shifts after interaction with the vortex. Section ^ shows how to derive the 
transverse force from the exact partial-wave solution of the Aharonov-Bohm 
problem for electrons. In our analysis the force on the vortex originates from 
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scattering of noninteracting quasiparticles. This is a valid assumption since 
normally the phonon mean-free path essentially exceeds the scale where the 
force arises (the phonon wavelength) EM . But in order to know the effect of 



this force on the whole superfluid, it is important to investigate how the force 
is transmitted to distances much larger than the mean-free path where the 
two-fluid hydrodynamics becomes valid. It is done in Sec. |[ The effect of 
the force at very large distances is also important for discussion of the Berry 
phase in Sec.^. The transverse force creates a circulation of the normal ve- 
locity at very large distances from an isolated vortex. Taking into account 
this circulation, the Berry phase yields a correct value of the transverse force, 
which agrees with the result derived from the momentum balance. 

The present paper does not discuss experimental aspects of this problem, 
which are addressed in other reviews devoted to rotating ^He and ^He P, ^ 

2 The Magnus force in classical hydrodynam- 
ics 

It is worth to remind how the Magnus force appears in classical hydrody- 
namics. Let us consider a cylinder immersed in an incompressible inviscous 
liquid. There is a potential circular flow around the cylinder with the velocity 

MO - ■ m 

Here r is the position vector in the plane xy, the axis z is the axis of the 
cylinder, and k is the circulation vector along the axis z. In classical hy- 
drodynamics the velocity circulation n = § ■ dl may have arbitrary values. 
There is also a fluid current past the cylinder with a transport velocity Vtr 
and the net velocity is 

v{r) = v^{r)+ Vtr ■ (2) 

This expression is valid at distances r much larger than the cylinder radius. 
At smaller distances one should take into account that the flow with the 
velocity vtr cannot penetrate into the cylinder, but velocity corrections due 
to this effect decrease as and are not essential for the further analysis. 
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The Euler equation for the hquid is 



^ + {v-V)v = --VP . (3) 
ot p 

Here p is the hquid density and P is the pressure. 

Assuming that the cyhnder moves with the constant velocity vl, i.e., 
replacing the position vector r by r — w^t, one obtains that 

— = -{vL-V)v. (4) 
Then the Euler equation yields the Bernoulli law for the pressure: 

^' = ^'o - ]j^p[v{r) - VlY = Pq- ^pvv{r)'^ - pv.,{r) ■ {vt^ - Vl) . (5) 

Here Pq and Pg = -Po~|p('^^tr — "^l)^ are constants, which are of no importance 
for the following derivation. Figure shows that due to superposition of two 
fluid motions given by Eq. (Q) the velocity above the cylinder is higher than 
below the cylinder. According to the Bernoulli law, the pressure is higher in 
the area where the velocity is lower. As a result of it, a liquid produces a 
force on the cylinder normal to the relative velocity of the liquid with respect 
to the cylinder. This is a lift force, or the Magnus force. 

In order to find the whole force, we must consider the momentum balance 
for a cylindrical region of a radius tq around the cylinder (see Fig. |ip. The 
momentum conservation law requires that the external force F on the cylin- 
der is equal to the momentum flux through the entire cylindrical boundary in 
the reference frame moving with the vortex velocity v^. The momentum-fiux 
tensor is 

Tlij = P5ij + pVi{r)Vj{r) , (6) 
or in the reference frame moving with the vortex velocity vl'- 

= P^ij + Pi^i - VLi){Vj - VLj) ■ (7) 

The momentum flux through the cylindrical surface of radius tq is given by 
the integral / dSjIl[j where dSj are the components of the vector dS directed 
along the outer normal to the boundary of the cylindrical region and equal 
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Figure 1: Magnus (lift) force. It is derived from the momentum balance 
inside a cylinder of radius vq. 
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to the elementary area of the boundary in magnitude. Then using Eqs. (|1|), 
(H), and the momentum balance yields the following relation: 

p[{vL - V,,) xk]=F . (8) 

On the left-hand side of this equation one can see the Magnus force as 
it comes in the classical hydrodynamics. The Magnus force balances the 
resultant external force F applied to the cylinder. In the absence of external 
forces the cylinder moves with the transport velocity of the liquid: vl = Vtr 
(Helmholtz's theorem). 

In the derivation we used the hydrodynamic equations only at large dis- 
tance from the vortex line, and the radius of the cylinder does not appear 
in the final result. Therefore Eq. (H) is valid even if the circular flow with 
circulation n occurs without any cylinder at all. In hydrodynamics such a 
flow pattern is called a vortex tube, a vortex line, or simply a vortex. Hydro- 
dynamics is invalid at small distance from the vortex line. This area is called 
the vortex core. But this does not invalidate the derivation of the Magnus 
force for a Galilean invariant liquid, in which the momentum is a well-defined 
conserved quantity even inside the vortex core where the hydrodynamic the- 
ory does not hold. 

In the momentum balance in a cylinder around a vortex, a half of the Mag- 
nus force is due to the Bernoulli contribution to the pressure, Eq. (^); another 
half is due to the convection term oc ViVj in the momentum flux. However, 
such decomposition of the resultant momentum flux onto the Bernoulli and 
the convection parts is not universal and depends on the shape of the area 
for which we consider the momentum balance. We may consider the mo- 
mentum balance in a stripe, which contain the vortex inside and is oriented 
normally to the transport flow (see Fig. ^). This yields again Eq. (|^), but 
now the pressure (the Bernoulli term) does not contribute to the transverse 
force at all, and only convection is responsible for the Magnus force. Then 
the physical origin of the Magnus force is the following. The liquid enter 
the stripe, which contains a vortex, with one value of the transverse velocity 
(equal to zero in Fig. |]) and exits from the stripe with another value of 
it. A difference between two values is Av in Fig. |^. The transverse force 
is the total variation of the transverse (with respect to the incident velocity 
v) liquid momentum per unit time. The latter is equal to a product of the 
current circulation p § Av ■ dl and the velocity v. 
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Figure 2: Momentum balance in a stripe area. The transverse force is deter- 
mined by the current circulation p § Av ■ dl round the stripe. 

3 The Magnus force in a superfluid 

In the superfluid state liquid motion is described by two- fluid hydrodynamics: 
the liquid consists of the superfluid and the normal component with the su- 
perfluid and the normal density Ps and and the superfluid and the normal 
velocity Vg and Vn respectively. The circular motion around the vortex line 
is related to the superfluid motion. Therefore Hall and Vinen P] suggested 
that the Magnus force is entirely connected with the superfluid density ps 
and the superfluid velocity Vs{r) = v^ir) + Vsn- Then instead of Eq. (||) one 
has: 

ps[{yL-Vs) ^ ti]= F . (9) 

Here and later on we omit the subscript "tr" replacing Vsu by Vs- But one 
should remember that the superfluid velocity Vs in the expression for the 
Magnus force is the superfluid velocity far from the vortex line. 

The "external" force F is in fact external not for the whole liquid, but 
only for its superfluid part. The force appears due to interaction with quasi- 
particles which constitute the normal part of the liquid, and therefore is 
proportional to the relative velocity vl — Vn 0- ^ Galilean invariant 
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liquid with axial symmetry the most general expression for the force F is 

F = -D{vl - Vn) - D'[z X {vl - . (10) 

The force component oc D' transverse to the velocity vl — Vn 'i^ possible be- 
cause of broken time invariance in presence of a vortex and resulting asym- 
metry of quasiparticle scattering by a vortex. 

Inserting the force F into Eq. (^), one can rewrite the equation of vortex 
motion collecting together the terms proportional to the velocity vl'- 

Pm[vl X k] + Dvl = p.s[Vs X k] + DVn + D'[z X Vn] . (11) 

The forces on the left-hand side of the equation are the effective Magnus force 
oc pm = ps — D'/k and the longitudinal friction force oc D. The forces on 
the right-hand side are driving forces produced by the superfluid and normal 
flows. In the theory of superconductivity the force Ps[vs x k], proportional 
to the superfluid velocity Vg, is called the Lorentz force. The left-hand side 
of Eq. (0) presents the response of the vortex to these driving forces. The 
factor Pm, which determines the amplitude of the effective Magnus force on 
the vortex, is not equal to the superfluid density ps in general. In the next 
sections we shall consider the contribution to D' from phonon scattering by a 
vortex (the lordanskii force). The contribution to D' from the bound states 
in the vortex core is discussed by Kopnin [|I4 . 



Nonlinear Schrodinger equation and two- 
fluid hydrodynamics 



In the Gross-Pitaevskii theory |^ the ground state and weakly excited states 
of a Bose g described by the condensate wave function ijj = aexp(z0) 

which is a solution of the nonlinear Schrodinger equation 

= -I^V^V^ + . (12) 

Here V is the amplitude of two-particle interaction. The nonlinear Schrodinger 
equation is the Euler-Lagrange equation for the Lagrangian 
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Noether's theorem yields the momentum conservation law 



where 



dt ^ ' 



dtp dt/j* 2 



(14) 
(15) 



is the mass current and 



dL 



dL 



2m 



(16) 



is the momentum-flux tensor. The pressure P in this expression corresponds 
to the general thermodynamic definition of the pressure via a functional 
derivation of the energy with respect to the particle density n = 



P = L 
where 



5E 

n- E 

on 



n 



dE 
dn 



V 



' dE ' 
dVn, 



-E 



|lV'l'-|;V=I^P, (17) 



E = -^^ H ^ip* 

dip dip* 



(18) 



is the energy density. But in the hydrodynamic limit (see below) the depen- 
dence of the energy on the density gradient is usually neglected. 

Using the Madelung transformation the nonlinear Schrodinger equa- 
tion (|l^) for a complex function may be transformed into two real equations 
for the liquid density p = ma? and the liquid velocity v = (fi;/27r)V0 where 
K = h/m is the circulation quantum. Far from the vortex line these equations 
are hydrodynamic equations for an ideal inviscous liquid: 



| + V(H=0. 



dv 
di 



(19) 
(20) 
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Here /i = Va?/m is the chemical potential. Equation ( ]T6| ) becomes the 
hydrodynamic momentum- flux tensor Iljj = P5ij + pVi{r)Vj{r). Thus the hy- 
drodynamics of an ideal liquid directly follows from the nonlinear Schrodinger 
equation. 

Suppose that a plane sound wave propagates in the liquid generating the 
phase variation 0(r, t) = (poexp^ik ■ f — icot). Then the liquid density and 
velocity are functions of the time t and the position vector r in the plane xy: 

p{f, t) = po + p(i) (f, t) , v{f, t) =Vo + {7(1) (f, t) , (21) 

where po and Vq are the average density and the average velocity in the liquid, 
P(i)(r, t) and 'U(i)(r, t) = are periodical variations of the density and 

the velocity due to the sound wave ((p(i)) = 0, ('^(i)) = 0). They should be 
determined from Eqs. (|19]) and (pOD after their linearization. In particular, 
Eq. (l20| ) gives the relation between the density variation and the phase 0: 

P(i) = ^m = -^^|^ + -o-V0(r-)| , (22) 

where Cs = ^JVaFJm is the sound velocity. Substitution of this expression 
into Eq. (p!9D yields the sound equation for a moving liquid with the wave 
spectrum uj = Cgk + k ■ Vq. 

The sound propagation is accompanied with the transport of mass. This 
is an effect of the second order with respect to the wave amplitude. The total 
mass current expanded with respect to the wave amplitude and averaged over 
time is 

J = PqVq + (P(l)'y(i)) + {p{2))Vo + po{Vi2)) , (23) 

where 

(P(i)^^(i)) = f\p) = Po<Pl^k = n{p)p . (24) 

is the mass current and n{p) is the number of phonons with the momentum 
p = hk and the energy E = e{p) + p ■ vq. The phonon mass current is the 
phonon momentum density in the reference frame moving with the average 
liquid velocity vo, and e{p) = CsP is the energy in the same reference frame. 

Mathematically the second-order corrections to the mass density, (p(2)), 
and the average velocity, {v(^2)), remain undefined, but there are physical 
constrains to specify them. First of all, we assume that phonon excitations 
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do not change the average mass density, and (p(2)) must vanish. As for the 
second-order correction {v(^2)) to the average velocity, it should produce a 
second-order correction (0(2)) to the phase which is impossible in quantum 
hydrodynamics. The simplest way to see it is to consider the propagation of 
phonons in an annular channel with the periodic boundary conditions. The 
phase variation over the channel length is a topological invariant, and weak 
excitations (phonons) cannot change it. Therefore ('^(2)) must vanish. More 
complicated arguments must be given for an open geometry, but intuitively it 
is clear that this basic physical constrain should not depend on the boundary 
condition. 

It is important to emphasize a difference between sound waves in a liquid 
and in an elastic solid. The sound wave in the elastic solid is not accom- 
panied by real mass transport: all atoms oscillate near their equilibrium 
positions in the crystal lattice, but they cannot move in average if the crys- 
tal is fixed at a laboratory table. Within our present formalism this means 
that the second-order contribution (^(2)) to the average velocity must not 
vanish, but compensate the second-order contribution (p(i)W(i)) to the mass 
current. Finally a phonon in a solid cannot have a real momentum but only 
a quasimomentum. The problem of the phonon momentum in liquids and 
solids has already been discussed a long time pOf (see also the recent paper 
by Stone |]30[), and they have noticed that the sound wave may have a dif- 



ferent momentum using Euler or Lagrange variables. In fact, the momentum 
should not depend on a choice of variables, but only on physical conditions. 
However, at various physical conditions a proper choice of variables can do an 
analysis more straightforward. In solids the Lagrange variables are preferable 
since in this case the velocity is related to a given particle and coincides with 
the center- of-mass velocity which must not change in average and therefore 
has no second-order corrections. Using Euler variables in liquids the aver- 
age velocity vq relates to a given point in the space and has no second-order 
corrections due to phonons. 

Expanding the momentum-flux tensor up to the terms of the second order 
with respect to the sound wave amplitude one obtains: 

Uij = Po6ij + poVoiVoj + U^l^ , (25) 

where the second-order phonon contribution is 

= {P(2))Sij + {p{l){V(^i))i)Voj + {P{l){V(^l))j)voi + Po{{V(i))i{V(^i))j) . (26) 
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The second-order contribution P(2) to the pressure can be obtained from the 
Gibbs-Duhem relation 5P = p6p at T = using expansions p = po + P(i) and 
p = Po + + /i(2), where /iq is the chemical potential without the sound 

2 

wave. This yields P(2) = PoP{2) + f^~2^! where |^ = Po/c^- According to the 
Euler equation @ the second-order contribution to the chemical potential is 

P{2) = — ^- Then 

(P \- n ^""W^ (27) 

(^(2))--^-po^. (27) 

In the presence of the phonon distribution the total mass current is 

/ = PqVq + j f^i.P) d3P = PoVo + ^ J n{p)pdsp , (28) 

In the thermal equilibrium at T > 0, the phonon numbers are given by the 
Planck distribution n{p) = no{E, Vn) with the drift velocity Vn of quasiparti- 
cles: 

^1 1 

no{E,Vn) = E{p)-p-v„ 7 ^ e{p)+p-{vo-Vn) 7 " (^^) 

exp ^^'j,^ — 1 exp ^ ' — — — 1 



Linearizing Eq. (|29| ) with respect to the relative velocity vq — Vn, one obtains 
from Eq. (pS]) that 

j = PqVq + Pn{Vn - Vq) . (30) 

This expression is equivalent to the two- fluid expression j = pvg + Pn{vn — 
Vs) = PsVs + PnVn assumiug that p = po = Ps + Pn, Vq = Vg, and the normal 
density is given by the usual two-fluid-hydrodynamics expression: 

1 f dnQ{e,0) 2 
Pn = -TTT^ / ^z. P dsP ■ (31) 



3/i3 J dE 

In the same manner one can derive the two-fluid momentum-flux tensor 0: 

Uij = PSij + PsVsiVsj + p (32) 

This analysis demonstrates that two-fluid-hydrodynamics relations can 
be derived from the hydrodynamics of an ideal liquid in presence of ther- 
mally excited sound waves, as was shown by Putterman and Roberts In 
order to obtain a complete system of equations of the two-fluid theory, one 
should take into account phonon-phonon interaction, which is essential for 
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the phonon distribution function being close to the equihbrium Planck dis- 
tribution. In the two-fluid theory the locally defined superfluid and normal 
velocities Vg and Vn correspond to the average velocity of a liquid in a fixed 
point of the space and to the drift velocity of the phonon gas respectively. 
The two-fluid hydrodynamics is valid only at scales exceeding the phonon 
mean-free path Ip^. 



5 Scattering of phonons by the vortex in hy- 
drodynamics 

Phonon scattering by a vortex line was studied beginning from the works by 
Pitaevskii and Fetter [RH] . We consider a sound wave propagating in the 



plane xy normal to a vortex line (the axis z). In the linearized hydrodynamic 
equations of the previous section the fluid velocity vq should be replaced by 
the velocity Vvir) around the vortex line: 



dt 



+ Po V ■ tT(i) 



Vp(i) 



dv, 



(1) 



dt 



+ f Vp(i) = - 
Po 



Using the vector identity 



(33) 
(34) 



for the velocity v = 
dt 



iv-V)v = V- 



Vy + {7(1), Eq. 



— V X [V X v\ 
41) can be rewritten as 



Po 



-V(t4 ■ %)) + [v(i) X K]S2ir) 



(35) 



(36) 



The perturbation from the vortex (the right-hand side) contains a 5-function 
because the vortex line is not at rest when the sound wave propagates 
past the vortex. In order to weaken the singularity one can introduce the 
time-dependent vortex velocity Vv{r,t) as a zero-order approximation for 
the velocity field 0. Then r in Eq. ([|) must be replaced by r — vit and 
dvv/dt = —{vl-'V)vv = —"V^Vl-Vv) + [vl x K]62{r)- Since there is no external 
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force on the liquid, the vortex moves with the velocity in the sound wave: 
vl = '^(i)(0,t). Now in the linearization procedure the fluid acceleration in 
Eq. ( PPD must be presented as dv/dt = dvy/dt + dv(i)/dt. As a result Eq. 
(|36D is replaced by: 



dv(^i) . c 



,2 



dt po 
Equation (|37|) yields 



+ ^Vp(i) = V[v, ■ {;(i)(f)] - ■ %)(0)] . (37) 



P(i) = -||^|^ + ^^.-[V0(r1-V0(O)]| . (38) 
Substitution of p(i) in Eq. ( p3D yields the linear sound equation for the phase: 



,,,-cfV^0--2.^.(f).V- 



(^1 - 



(39) 



In the long-wavelength limit k one may use the Born approximation. 
The Born perturbation parameter nk/cg is on the order of the ratio of the 
vortex core radius ~ n/cg to the wavelength 27r/fc. Then after substituting 
the plane wave into the right-hand side of Eq. ( |5UD the solution of this 
equation is 



(f) = exp(—iiut) |exp(iA; ■ r) 



"i7 / d2rM^\k\f-n\)k-v,{f,)[2exp{ik-f,) - 1]| . (40) 

Here H^^^ (z) is the zero-order Hankel function of the first kind, and ^HQ^\k\f— 
ri|) is the Green function for the 2D wave equation, i. e., satisfies to the 
equation 

-(A;2 + V2)0(r) = 52(r-fi) . (41) 

In the standard scattering theory they use the asymptotic expression for 
the Hankel function at large values of the argument: 

lim hI.^\z) = J Ae^(--/4) . (42) 
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If the perturbation is confined to a vicinity of the hne, then ri ^ r and 



\f—fi\ ^ r . (43) 



After integration in Eq. (^Of ) the wave at fcr 3> 1 becomes a superposition of 
the incident plane wave oc exp{ik ■ r) and the scattered wave oc exp(i/cr): 



= 0oexp(— 



exp[tk ■ r) H — exp[tKr) 



(44) 



Here a{(p) is the scattering amphtude which is a function of the angle ip be- 
tween the initial wave vector k and the wave vector k' = kf/ r after scattering 
(see Fig. |^). In the Born approximation (see the paper and references 
therein for more details) 

, . [k 1 ji.^ -1 / \ 1 jk K .^r sini/Jcosv? 

a((p) = \ e* 4 k X k ] ■ k— 1 - -Vtt = — \ * , (45) 

V27rc, ^ 2k^J 2 V 27r c, l-cosy? ' ^ ^ 

where q = k — k' is the momentum transferred by the scattered phonon to 
the vortex, and = 2k^{l — cosy?). 

Thus the vortex is a line defect which scatters a sound wave. Scattering 
produces a force on the defect (vortex). If the perturbation by the line defect 
is confined to a finite vicinity of the line, the force 

FP'' = aiicsf'' - a±Cs[z x jP''] (46) 

is determined by two effective cross-sections P, 0]: the transport cross-section 
for the dissipative force, 

(T|| = J a{(p){l — cos(p)dip , (47) 
and the transverse cross- section for the transverse force, 

a± = a {(p) sin (fdif . (48) 



The differential cross-section a^ip) = |a(v9)p with a{(p) from Eq. (|45|) is 
quadratic in the circulation k, and even in ip. Therefore in the Born approx- 
imation the transverse cross-section a± vanishes. 
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Figure 3: Scattering of a sound wave by a vortex. 



However, due to a very slow decrease of the velocity oc 1/r far from 
the vortex, the scattering amplitude is divergent at small scattering angles 
^ 0: 

lim a(o9) = -W e *—. (49) 

This divergence is integrable in the integral for the transport cross-section, 
Eq. P7|), and its calculation is reliable. Contrary to it, the integrand in 
Eq. (^81) for the transverse cross-section has a pole at </9 = 0. The principal 
value of the integral vanishes, but there is no justification for the choice 
of the principal value, and the contribution of the small angles requires an 
additional analysis. 

At small scattering angles <^ l/\/kr the asymptotic expansion Eq. 
( ^31) is invalid, and one cannot use the differential cross-section or the scat- 
tering amplitude for description of the small-angle scattering 0, |^, |, ^ 



Meanwhile, the small-angle behavior is crucial for the transverse force as 
demonstrated below. The accurate calculation of the integral in Eq. (EOf) for 
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small angles (see Ref. and Appendix B in Ref. ^) yields that at ^ 1 



I exp{—iujt + ik ■ r) 



Using an asymptotic expression for the error integral 



dt — ^ — 



exp(— 2; 



(50) 



(51) 



at z 



00, one obtains for angles 1 ^ (f ^ l/y/kr: 

ink 



I exp(— 2u;t) 



exp(iA; ■ r ) 1 + 



IK 



k 1 



2cs \^\ 



Cs V 27rr if 



exp ( ikr + 



(52) 

The second term in square brackets coincides with scattering wave at small 
angles -C 1 with the amplitude given by Eq. (|49|). But now one can 
see that the standard scattering theory misses to reveal an important non- 
analytical correction to the incident plane wave, which changes a sign when 
the scattering angle crosses zero. Its physical meaning is discussed in the 
next section. 



6 The lordanskii force and the Aharonov— 
Bohm effect 

The analogy between the phonon scattering by a vortex and the Aharonov- 
Bohm effect for electrons scattered by a magnetic-flux tube becomes evident 
if one rewrites the sound equation (|39| ) in presence of the vortex as 




fcV - -^V + —vA = 0. (53) 



It differs from the sound equation (|39D by the term of the second order 



m. Vy (X K and by absence of the contribution from the vortex-line motion 
[the term oc 0(0) on the right-hand side of Eq. (^)]. This difference is 
unimportant for the calculation of the transverse force, which is linear in k. 
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On the other hand, the stationary Schrodinger equation for an electron in 
presence of the magnetic flux $ confined to a thin tube (the Aharonov-Bohm 
effect §) is: 



1 

2m 



-ihV - -A] il){r) 



(54) 



Here if) is the electron wave function with energy E and the electromagnetic 
vector potential is connected with the magnetic flux $ by the relation similar 
to that for the velocity around the vortex line [Eq. (Hll 



A = $ 



X r] 
27rr2 



(55) 



Let us consider the quasiclassical solution of the sound equation: 



()o exp —iujt + ik ■ r + 



I exp{—iujt + ik ■ f) ( 1 + 



ink 
2'7rc, 



(56) 

where SS = —{hk/cg) v^-dl = h9nk/2TiCs is the variation of the action due 
to interaction with the circular velocity from the vortex along quasiclassical 
trajectories. The angle 9 is an azimuth angle for the position vector r mea- 
sured from the direction opposite to the wave vector k (see Fig. |). This 
choice provides that the quasiclassical correction vanishes for the incident 
wave far from the vortex. One can check directly that Eq. (|56|) satisfies 
the sound equation (^) in the first order of the parameter nk/cg. For the 
Aharonov-Bohm effect the phase 5S/h arises from the electromagnetic vector 
potential: 5S = — (e/c) f A - dl. 

The velocity generated by the sound wave around the vortex is 



V{i) = — V0 = —00 exp(— icut + ik ■ f) [ik 



27r 



2tt 



(57) 



From Eqs. (|38|) and (^7\) we can obtain the phonon mass current with first- 
order corrections in v^. 



Jph = {P(1)V{1)) 



1 



dkn{hk)h 



k k 
k v^ir) - {k ■ v^{r)) — 

Cq Cq rv 



(58) 



Due to the last term in this expression, the phonon mass current is not 
curl-free. 
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5S+ 




Figure 4: Aharonov-Bohm interference. 

According to Eq. (p6D the phase is multivalued, and one must choose a 
cut for an angle 6 at the direction k, where 6 = ±7r . The jump of the phase 
on the cut line behind the vortex is a manifestation of the Aharonov-Bohm 
effect 1^: the sound wave after its interaction with the vortex has different 
phases on the left and on the right of the vortex line. This results in an 
interference 0. 

In the interference region the quasiclassical solution is invalid and must 
be replaced by Eq. (|50D . The width of the interference region is dint ^ 
\Jr /k. Here r is the distance from the vortex line. The interference region 

corresponds to very small scattering angles ~ dint/r = l/\/kr. 

Now we are ready to consider the momentum balance using the condition 
that / dSjll^j = where subscript _L is for a component normal to the 
wave vector k of the incident wave. The total momentum-flux tensor can be 
obtained from Eqs. (^) and (|26|) assuming Vo{r) = v^ij) +Vs and neglecting 
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some unimportant terms: 

Hjj = -po{Vs - Vl) ■ Vjij + PoVoiVoj 

+ {P{i)iv(i))i)v^j + (P(i)(t;(i))j>« + Po((^^(i))i(^^(i))i) • (59) 

The first two terms in this expression yield the momentum flux without 
phonons, which gives the Magnus force for a hquid with the density po and 
the velocity Vg. The term oc v^j does not contribute to the momentum 
flux through a cylindrical surface around the vortex, since the velocity is 
tangent to this surface. The term Vy.{p[i){v(^i))j) , in which the mass current 
is given Eq. ( p4|) for the plane wave in absence of the vortex, gives 
a contribution to the momentum flux, which exactly cancels the contribution 
of the term Po(('^(i))i('^(i))i) outside the interference region, where the velocity 
is given by Eq. (^7|). Finally only the interference region contributes the 
momentum flux / dSjIl±j. 

In the interference region the velocity is obtained by taking the gradient 
of the phase given by Eq. (|50| ) . Its component normal to k, 

V(i)± = - — — = (poexp{-iujt + ik-r)- ^— , (60) 

Zirr oip AttCsT 0(p 

determines the interference contribution to the transverse force: 
J dSjPo{iv(^i))±iv(^i))j) = J po((f(i))±(f(i))r)rc?¥P 
^^Mt^lSS_-6S4^J^p,K^^^.f\ (61) 

where 6S± are the action variations at ^ =i=7r. Thus the interference region, 
which corresponds to an inflnitesimally small angle interval, yields a flnite 
contribution to the transverse force, which one could not obtain from the 
standard scattering theory using the differential cross- sect ion. In fact the 
details of the solution in the interference region are not essential: only a 
jump of the phase across the interference region is of importance. 

Finally the momentum balance condition / dSjIl±j = yields the relation 

PolivL - Vs) X /€] - [f\p) xk]=0. (62) 
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The second vector product on the left-hand side is a transverse force which 
corresponds to the transverse cross-section [see Eq. (|^)] 



5S- - 5S+ 

equal to n/cg. For the Planck distribution of phonons, j^^ must be replaced 

by pni^n - Vs)- 

Ps[{vl - Vs) X k]+ Pn[{vL - Vn) X = . (64) 

The term oc [vi — Vn) is the lordanskii force, which corresponds to D' = —npn 
in Eq. (0). The longitudinal force oc D is not present in Eq. (|^) since we 
ignored terms of the second order in k in order to simplify our derivation. 
According to Eq. (|^) the vortex moves with the center-of-mass velocity 

V = ^Vs + ^Vn- 

p A I p II- 

Our scattering analysis was done in the coordinate frame where the vortex 
is at rest and we neglected the relative superfluid velocity Vg — vi with respect 
to the vortex. If the velocity Vs — vl is high it can affect the value of the 
effective cross-section. But since the phonon momentum is linear with respect 
to Vn — Vs, the dependence of the cross-section on Vs — vl is a. nonlinear 
effect. Thus our momentum balance took into account all effects linear in 
the superfluid velocity Vs- This is confirmed by a more elaborate analysis of 
Stone m . 



7 Partial-wave analysis and the Aharonov— 
Bohm effect 

Studying interaction of phonons with a vortex we solved the sound equa- 
tion using the perturbation theory. It is completely justified because the 
parameter of the perturbation theory nk/cs is the ratio between the vortex 
core and the phonon wavelength, which is always small for phonons. But in 
the Aharonov-Bohm problem for electrons the corresponding parameter is 
7 = $/$i, where $i = hc/e is the magnetic-flux quantum for one electron 
(two times larger than the magnetic-flux quantum $o = hc/2e for a Cooper 
pair). This parameter can be arbitrary large, and the perturbation theory is 
not enough to describe an expected periodic dependence on 7. On the other 
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hand, there is an exact solution of the Aharonov-Bohm problem for electrons 
obtained by the partial- wave expansion 0, and it will be shown now how 
to derive the transverse force from this solution. Another derivation of the 
force on the Aharonov-Bohm flux tube using the wave-packet presentation 
is discussed by Shelankov [|12| . 

We look for a solution of Eq. (|5^) as a superposition of the partial 
cylindrical waves = X) '^/(^) exp(iZy9). Partial-wave amplitudes ipi should 

satisfy equations in the cylindrical system of coordinates {r,ip): 

-TT J 5 — + k ipi = Q . (65) 

Here k is the wave number of the electron far from the vortex so that E = 
h'^k'^/2m. We need a solution of Eq. (pSI), which at large distances has an 
asymptotic behavior given by Eq. 



2n 

Vi = y exp 



.vr , _ 
2 



TT 71 

cos{kr - -I + 5i - -) , (66) 



nkr 

where n is the particle density and the partial-wave phase shifts are 

<5^ = (/-|/-7|)7r/2. (67) 



For Si = Eq. ( pB| ) yields the partial-wave amplitudes of the incident plane 
wave at large distances r. But for nonzero 6i there is also the scattered wave 
in Eq. (H3) with the scattering amplitude 



" V i Gl) ~ exp{2i6i)] exp{ilip) . (68) 
The transverse force is determined by the transverse cross-section 

(7_L = J \a((f)\'^ sm(fd(f 



Shifting the number / by 2 in the first sum and by one in the fourth sum 
one obtains the expression for the transverse cross-section in the partial-wave 
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method derived long ago by Cleary [^: 



cr± = y" |«(V5) P sin Lfdi^ sin 2(5; - 5i+i) . (70) 



Using the phase shift values for the Aharonov-Bohm solution, Eq. (p7|), the 
transverse cross-section is 

(T_L = — — sin27r7 . (71) 

K 



However, the shift of / in the first sum of Eq. ( pQ]) is not an innocent operation 
because of divergence of the first and second sum at / — ±oo. The derivation 
of Clearly' s formula ( [70| ) assumes that the first and the second divergent sums 
in Eq. (^) should exactly cancel after the shift of /. But if one does not 
shift /, the difference of the first and the second sum is finite. Moreover, this 
difference cancels the contribution of the third and the fourth sum, and a± 
vanishes in the first order with respect to the phase shifts 5i. 

Ambiguity in the calculation of the partial-wave sum for the transverse 
cross-section is another manifestation of the small-scattering-angle problem 
in the configurational space: the standard scattering theory does not provide 
a recipe how to treat a singularity at zero scattering angle. The zero- angle 
singularity is responsible for a divergent partial-wave series. The way to avoid 
the ambiguity is similar to that in the configurational space: one should not 
use the concept of the scattering amplitude for calculation of the transverse 
force. 

We must analyze the momentum balance. The momentum-flux tensor for 
the electron Schrodinger equation (Q) is 



-'f (^-ihV, - -^A^ (-^^Vj- - -^A^ Vj-t/'I . (72) 

If the axis x is directed along the wave vector k of the incident wave, 
the transverse force is determined by the momentum-fiux tensor component 
Ilyri where r is the radial coordinate in the cylindrical coordinate system 
r = ^/x"^ + y"^ and if = arctan(?//x). In Eq. ([T^) the terms A and d/rd(p, 
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which are inversely proportional to r, are not important at large r. Then 



n„r = Re < sin 09 



e*('-'>[> , (73) 



dr dr dr 
where ipi are given by Eq. (|66|). Finally the transverse force is 




(74) 



Here we also made a shift of / in some sums, but it is not dangerous because 
the terms of these sums decrease in the limits / ±oo. 

Inserting the partial- wave amplitudes and their radial derivatives into the 
expression for the force we obtain: 

Fj_ = -kn— V sin 2((5i+i - 5i) = -jva^ , (75) 
m ^ 

where j = hkn and v = hk/m are the momentum density and the velocity in 
the incident plane wave, and a± is the effective cross-section Eq. (^) derived 



by Cleary ||36[ . Thus there is a well-defined effective transverse cross-section, 
despite that we cannot use directly its expression (^8|) via the differential 
cross-section. But formally we can "repair" this expression by a recipe how 
to treat the singularity at small angles: we should add to the differential 
cross-section a singular term proportional to the derivative of S{ip), and to 
take the principle value of the integral over the rest part of the differential 
cross-section. 

One can obtain the transverse cross-section for phonons directly from Eq. 
( |70D assuming that sin27r7 ^ 27ij = —Kk/cg. However, it is useful to follow 
the connection between Cleary's formula (|70D and Eq. (^) obtained from the 
quasiclassical solution. In the classical limit the partial wave I corresponds 
to the quasiclassical trajectory with the impact parameter b = l/k = M/p 
and the small scattering angle is = —d6i/dl = —{l/2hk)d6S{b)/db, where 
6S{b) is the action variation along the trajectory with impact parameter b. 
Finally replacing the sum by an integral Eq. (|70| ) can be rewritten as 



k J-oo dl hk J-oo db 
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which yields (^) since 6S±. = ^^(±00). Strictly speaking Eq. ( |76D is valid if 
6S{b) is a continuous function. This is the case for rotons which are scattered 
quasiclassically |§, 01- But it yields a correct answer even for phonons de- 
spite 6S{b) has a jump at 6 ~ 0. Thus even though phonon scattering cannot 
be described in the quasiclassical approximation (it yields the scattering an- 
gle ip = —{l/hk)d5S{h)/dh ~ at 6 7^ 0), the quasiclassical expression (ff^) 
gives a correct phonon transverse cross-section. 



8 Momentum balance in the two-fluid hydro- 
dynamics 

Up to now we analyzed spatial scales much less than the mean-free-path 
Iph of phonons (ballistic region). Sound waves (phonons) interacted with 
the velocity field generated by a vortex, but phonon-phonon interaction was 
neglected. Now we shall see what is going on at scales much larger than Iph, 
where the two-fluid hydrodynamics is valid. 

Scattering of phonons by the vortex in the ballistic region produced a 
force F on the vortex, and according to the third Newton law a force —F 
on the phonons (the normal fluid) should also arise. A momentum transfer 
from a phonon to a vortex takes place at distances of the order of the phonon 
wavelength A = 27r/fc, which is much less than the mean-free path Iph- This 
means that at hydrodynamic scales the force on the normal fluid is a 5- 
function force concentrated along the vortex line. The response of the normal 
fluid to this force is described by the Navier-Stokes equation with the dynamic 
viscosity ?7„: 

^TT + (^n ■ V)fn = Vnl^Vn VT , (77) 

dt p PnP 

where z/„ = rjn/ Pn is the kinematic viscosity, S is the entropy per unit volume, 
and T is the temperature. Our problem is similar to the Stokes problem of a 
cylinder moving through a viscous liquid [0] . Neglecting the nonlinear inertial 
(convection) term {vn ■ V)vn in the Navier-Stokes equation, the (5-function 
force produces a divergent logarithmic velocity field (the Stokes paradox ||l|]): 

F r 

Vn{r) =Vn + In -r- , (78) 

47r?7„ Iph 
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where Vn is the normal velocity at a distance of order r ^ Iph, which separates 
the ballistic and the hydrodynamic regions. In this expression the mean-free 
path Iph replaces the cylinder radius of the classical Stokes problem in the 
argument of the logarithm. Such a choice of the lower cut-off assumes that 
the quasiparticle flux on the vortex is entirely determined by the equilibrium 
Planck distribution at the border between the ballistic and the two-fluid- 



hydrodynamics region at r ~ Iph [|37 



However small the relative normal velocity Vn — vl could be, at distance 
of the order or larger than Vm ^ i^n/lvn — Vl\ the nonlinear convection term 
becomes important and stops a logarithmic growth of the normal velocity. 
Due to the force F the normal velocities Vnoo and Vn at large (r ~ Vm) and 
small (r ~ Iph) distances from the vortex line are different (the viscous drag 

i): 

F= , i^'^Jj A vnoo-v^). (79) 
ln(r^//p;,) 

At very large distances r ^ Tm the nonlinear convection term is more impor- 
tant than the viscous term. Thus the scale r^, which we shall call Oseen's 
length, divides the hydrodynamic region onto the viscous and convection 
subregions [^. All relevant scales are shown in Fig |[ For a longitudinal 



force the solution of the Navier-Stokes equation valid for both the viscous 
and convection subregions was obtained by Oseen long ago (see Ref. [IH). 

The force —F is transmitted to large distances by a constant momentum 
flux. In the viscous subregion momentum transport in the normal fluid occurs 
via viscosity: Fi = — § TijdSj, where Tij = —rjn (ViVnj + ^ jVni) is the viscous 
stress tensor. On the other hand, the total momentum flux for the whole 
liquid should vanish: / liijdSj + / TijdSj = 0, where Uij is given by Eq. (|3^). 
The normal velocity field does not contain the circular velocity Vy, and there 
is no normal circulation in the viscous subregion. Therefore the momentum 
flux fUijdSj yields the superfluid Magnus force, i.e., the momentum balance 
in the viscous subregion confirms that the force F satisfies Eq. (^. 

In the convection region the viscosity becomes ineffective and the momen- 
tum conservation gives again fUijdSj = 0, like in the ballistic region. The 
superfluid part of the momentum flux is related to the superfluid Magnus 
force oc ps and the normal part of the flux should transmit the same force F. 
The relation between the force F and the relative normal velocity Vnoo ~ vl 
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Figure 5: Relevant scales: the phonon wavelength A, the phonon mean free 
path Iph, and Oseen's length ~ T^n/\vn — vl\ 
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in the convection subregion can be derived from Eqs. ( jlOD and (fTOl): 



F = -D{vl - Vnoo) - D'[Z X {vl - Vnoo)] , (80) 

where the parameters D and D' are connected with D and D' in Eq. (p!0|) 
by a complex relation 

1 ln(r„//pft) 1 



D + iD' 4:7Tr]n D + iD' 

Since in the convection region the viscosity is ineffective and the normal fluid 
behaves as an ideal one, the only way to transmit the transverse component 
of the force F to infinity is to create a circulation of the normal velocity P^]. 
Because of the viscous drag separation on the longitudinal and the transverse 
force should be done with respect to the normal velocity Vnoo — vl, but not 
Vn — vl- Thus the normal circulation is determined by D', but not D': 

dl- v„, = -— = z o , (82) 



Pn 1 _^ 



npn ln(rm/lph) 



where we neglected the longitudinal component cx: Z) in the ballistic region 
and used the relation D' = —npn for the lordanskii force. 

In the convection subregion more problematic is to transmit not the trans- 
verse, but the longitudinal component of the force. A longitudinal force is 
impossible if viscosity is neglected completely: a body moving through an 
ideal liquid does not produce any dissipative force (d'Alembert's paradox 
[^). The paradox is resolved by finding that a laminar wake should arise 
behind a moving body. Within the wake one cannot neglect viscosity even 
deeply in the convection area when r 3> r^. The width of the laminar wake is 
growing as ~ i/rr^ far from the moving body [Q. Solving the Navier-Stokes 
equation by Oseen's method one can find out how the laminar wake and the 
normal circulation are formed during the crossover from the viscous to the 
convection subregion. One can find a detailed analysis of this crossover in 



the presence of the longitudinal and the transverse force in Ref. ^ . 

The force F is the mutual friction force introduced by Hall and Vinen 
for the analysis of propagation of the second sound in rotating superfiuid. 
With help of Eqs. (|^), (|10D , and (|79| ) one obtains a linear relation between 
the mutual friction force and the counterfiow velocity Vs — Vnoo bearing in 
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mind that the average normal velocity practically coincides with the velocity 

F = B[Vs - Vnoo) 2p~ ~ ' ^ ' 

The Hall-Vinen parameters B and B' are given by a complex relation [^] 

2p 1 1 ^ ^ 1 ^ 1 ^ Hrjlph) ^ 1 .g^. 

KpsPn B - iB' ipsK D + iD' ips^, 47r?7„ D + iD' ' 

According to this relation a strong viscous drag (large logarithm ln(rm/0) 
or small viscosity rj^) suppresses the effect of the transverse force (oc D') on 
the mutual friction [^. But the effect of the superfiuid Magnus force and the 
longitudinal force oc D is also suppressed in this limit. In fact, this is a limit 
when the force from quasiparticle scattering (transverse, or longitudinal) is 
so strong that the normal fluid sticks to the vortex, and Vl = Vn- Then 
the resultant force F depends only on viscosity, as in the classical Stokes 
problem. 



9 Magnus force and the Berry phase 

Now let us consider a connection between the transverse force and the Berry 
phase. We shall use the hydrodynamic description with the Lagrangian ob- 
tained from Eq. (|13|) after the Madelung transformation: 

L = - ^V0^ - ^p^ . (85) 

The first term with the first time derivative of the phase (p (Wess-Zumino 
term) is responsible for the Berry phase G = ASb/^i, which is the variation of 
the phase of the quantum-mechanical wave function for an adiabatic motion 
of the vortex around a closed loop |]l6l . Here 



ASb = I dfdt = - 1 dfdt ^{VL ■ VlH . (86) 

is the classical action variation around the loop and Vl4> is the gradient of the 
phase 0[r — rL(t)] with respect to the vortex position vector rL(t). However, 
Vl0 = — V0, where V0 is the gradient with respect to r. Then the loop 
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integral § dl yields the circulation of total current j = {h/2tx) < pV(j) > 
for points inside the loop, but vanishes for points outside. As a result, the 
Berry-phase action is given by [0 

ASB = V^fidl-j) . (87) 

where V is the volume inside the loop (a product of the loop area and the 
liquid height along a vortex). Contrary to Eq. (|86D , the loop integral in Eq. 
( P?!) is related with the variation of the position vector f, the vortex position 
vector being fixed. 

Since the Berry phase is proportional to the current circulation, which 
determines the transverse force, there is a direct connection between the 
Berry phase and the amplitude of the transverse force on a vortex, as was 
shown, e.g., in Ref. |T^.Thus in order to find the Berry phase, one should 
calculate the current circulation. If the circulation of the normal velocity 
at large distances vanished (as assumed in Refs. |^ 0, |TB|), the current 
circulation would be f{dl ■ j) = psK,, and the Berry phase (as well as the 
effective Magnus force) would be proportional to ps. However, according to 



Sec. H and and the recent analysis by Thouless et al. at very large 

distances (in the convection subregion) the normal circulation k„ = §{dl -Vn) 
does not vanish and §{dl ■ j) = psK + pn^n- Using a proper value of the 
asymptotic normal circulation Kn given by Eq. (|8^) the Berry phase yields 
the same transverse force as determined from the momentum balance. 

In order to obtain a correct value of the transverse force from the Berry 
phase, one should choose a loop radius much larger than Oseen's length r^- 
If the loop radius is chosen in the viscous subregion Ip^ < r < r^, the total- 
current circulation is proportional ps, but the Berry phase does not yield 
the total transverse force, since a part of it, namely, the lordanskii force, is 
presented by the viscous momentum flux, which cannot be obtained in the 
Lagrange formalism. And if the loop radius is chosen in the ballistic region, 
the total-current circulation is not defined at all and depends on a shape of 
the loop, since the phonon mass current is not curl-free, as pointed out after 
Eq. m. 
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10 Discussion and conclusions 



The momentum-balance analysis definitely confirms an existence of the trans- 
verse force on a vortex from phonon scattering (the lordanskii force). This 
conclusion agrees with the results of the recent analysis of Thouless et al. 

The Berry phase yields the same value of the transverse force as the mo- 
mentum balance, if a proper value of the normal circulation at large distances 
from the vortex is used for the calculation of the Berry phase. However, the 
Berry-phase analysis itself cannot provide the normal-circulation value, since 
the latter is determined by the processes at small distances from the vortex, 
which are beyond of the Berry-phase analysis. The small- distance processes 
determine the force between the superfluid and the normal component, which 
is present in the small-distance boundary condition for the Navier-Stokes 
equation in the two-fluid-hydrodynamics region. The transverse force on the 
normal component is transmitted to infinite distances by the constant mo- 
mentum flux, which requires a normal circulation far from the vortex. This 
circulation should be used for determination of the Berry phase. 

Ambiguity of the Berry-phase analysis of the transverse force originates 
from ambiguity of the transverse force in the Lagrange formalism, which was 
discussed in the end of Sec. V in Ref. 0). Adding a constant C to the 
total density in the Wess-Zumino term, i.e., replacing p by p + C, one obtain 
a different amplitudes of the Berry phase and the transverse force without 
any effect on the field equations for the condensate wave function (see also 
discussion of this constant in Fermi liquids by Volovik This arbitrari- 

ness has a profound physical meaning. Derivation of the Magnus force in the 
Lagrange formalism dealt only with large distances much exceeding the vor- 
tex core size. However, the processes inside the core affect the total Magnus 
force in general. Deriving the Magnus force from the momentum balance we 
use the condition that the total momentum flux / IlijdSj through the sur- 
face of large radius around the vortex line vanishes. This is true only for a 
Galilean invariant liquid satisfying the momentum conservation law. If the 
liquid in the vortex core interacts with the external world (e.g., with crystal 
impurities in superconductors), the momentum balance condition must be 
§ liijdSj = fi where / is the force on the vortex core which could also have 
a transverse component (the Kopnin-Kravtsov force [|14|)- One can calculate 
such a force only from the vortex-core analysis. The results of this analysis 
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must be used for determination of an unknown constant in the Wess-Zumino 
term. 
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